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A complete characterization of quantum fluctuations in many-body systems is accessible through 
the full counting statistics. We present an exact computation of statistical properties of light in a 
basic model of light-matter interaction: a multimode photonic field coupled to a single two-level 
emitter. We mostly consider an initial coherent state in a given mode and demonstrate how the 
original Poissonian statistics gets modified because of quantum many-body scattering effects leading 
to non-Poissonian distributions. We argue that measuring this statistics in a simple quantum optical 
setup provides an insight into many-body correlation effects with photons. 



Introduction. — The concept of the Full Counting 
Statistics (FCS) has been introduced first in the 
quantum-optical context in order to character- 

ize statistical properties of a non-interacting quantized 
electromagnetic field. Its fully quantum derivation has 
been presented in Ref. [3]. Later on, this concept 
has been borrowed and actively developed in the field 
of mesoscopic physics [4 for studying statistical prop- 
erties of electronic currents in meso- and nanoscopic 
devices for non-interacting and interacting electrons 
[5],I5],[7],[g],[i],[TD],[II|. It has been recently shown 
that the FCS has relationships to a classical-quantum 
crossover description [12], a quantum entanglement [13], 
a characterization of phase transitions [2]. The FCS of 
nonlocal observables can be used to quantify correlations 
[15H17j and a prethermalization behavior in many-body 
systems [H], as well as to define a certain topological 
order parameter |19j . 

Motivated by these developments we revisit the origi- 
nal problem of computing the FCS for photons interact- 
ing with an emitter. We give it the full quantum consid- 
eration treating the interaction nonperturbativcly. This 
study on the interface of quantum optics and nanoscopic 
physics becomes very actual nowadays, when the two 
fields are merging together in the continuous progress of 
fabrication of low- dimensional hybrid photon-solid state 
nanodevices [20]. One of the objectives in this interdis- 
ciplinary research is to obtain strong photon nonlinear- 
ities as well as a strong photon-emitter interaction for 
the purposes of an efficient control over individual atoms 
and phonons. In this respect, the knowledge of statis- 
tical properties of an interacting photon-emitter device 
becomes essential. 

In this Letter we present an exact calculation of the 
FCS in a basic model of light-matter interaction (see 
Fig.jlj); a multimode propagating photonic field interact- 
ing with a two-level emitter. A number of recent studies 
have focused on a low-dimensional transport properties 
of a few-photon initial states [H]. Here we significantly 
generalize and extend these studies in two respects. First 
we consider rather generic initial states (while mainly fo- 
cusing on experimentally-easy realizable coherent state). 
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FIG. 1: Our system consists of a two-level emitter situated 
at a; = and coupled to a waveguide (transmission line). 
A pulse of the length L is prepared in the initial coherent 
state lofco) in the mode fco. We assume the linear spectrum 
for the radiation field. Because of the interaction with an 
emitter photons can scatter both forward and backward. The 
distribution of photons after scattering can be detected by 
means of existing experimental techniques. 



Second, we assume this initial state prepared as a pulse 
of a size L which introduces an important realistic pa- 
rameter in photon's statistics. Third, going beyond first 
and second-order coherences we compute the whole FCS 
which is a source of all quantum statistical correlations 
hidden in this fundamental interacting system. Many- 
body correlations are present here because of a mul- 
timode structure of the photonic field interacting with 
emitter. Interactions considerably modify the statistics 
of photons in forward and backward scattering channels 
in comparison with the Poissonian one of the incident co- 
herent light. In a particular limit of continuous laser field 
(large L) we observe several universal features, e.g. for- 
mation of bimodal distribution for arbitrary initial state. 
We show how tuning the parameters of initial state and 
pulse duration reveals many-body effects. This can be 
useful for control of and manipulations with quantum 
states in emerging nanophotonic devices and circuits. 

The system. — The most basic model of light- matter in- 
teraction is described by the Hamiltonian H — Hph + 
Hem + -Hp/i-em, whcrc the free photon part, Hph — 
I f^fc C ajfcO,{fej features the linear, unbounded dis- 
persion with two branches ^ = ± = r, / (right and 
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left); the Hamiltonian of an emitter is approximated 
by a two- level system with the transition frequency A, 
Hem = ycr^; and the photon-emitter interaction is 
treated in the dipole and the rotating wave approxima- 
tions, Hph-em = 9oJ2^I dk{a'^^a~ + a^kcr~^). As usual, 
cr* and (7^ denote the Pauli matrices. Note that we 
consider the one-dimensional waveguide geometry, and 
therefore the transverse mode index is suppressed. 

Introducing symmetric Uek = (arfe + a/,-fc)/\/2 and 
antisymmetric Uok = [o-rk — o-i-k)I^P^ held configura- 
tions, we decouple the initial Hamiltonian into the sum 
H = Hp -|- iJ„ of the ewen, 



even-odd basis 
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dk 



k a 



k^-ek + g 
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where g — go\/2, and the orfd. Ho = J dk k alf,aok, con- 
tributions. For a local scatterer, one can always choose 
such a basis in which only one, even, mode is scattered, 
while the other, odd, mode is decoupled from the in- 
teraction Hamiltonian and therefore propagates without 
scattering. Both even and odd photons are additionally 
labelled by momentum values k lying on a single branch 
of the linear dispersion (see the Supplementary Material 
[22] for details). 

The major consequence of the even-odd decoupling is 
a factorization of the scattering matrix S into the prod- 
uct S = SeSo, where So equals the identity operator. 
The matrix Se corresponding to the model ([T]) has been 
constructed in Ref. |23j for all sectors of photon numbers. 
This construction allows for the full quantum description 
of scattering of an arbitrary initial field configuration in- 
cluding the coherent light. The scattering of the latter 
has been also studied in Ref. [25] . 

In order to completely specify the model presented in 
Fig.jlj we define the initial state to be |in) = |Q;)r,fco'8)| i), 
where the incident right-moving photons are prepared 
in the coherent state \a)r.ko = \ct)r = -Dr,fco(Q^)|0), and 
the two-level emitter is initially in the ground state | l) . 
Here |0) denotes the photonic vacuum, and Dr.ka{ct) = 
Dr{a) — exp(Q;6]!j,^ — Q;*6r.feo) is the displacement oper- 
ator. The latter is composed of the operators b\ = 

bt = ^ '^2;aj!(x)e*'^°'^ creating wave packets cen- 

tered around the given mode ko of the right branch of 
the original spectrum with the width 2'k/L (see the Sup- 
plementary Material [53] for details). Analogously we 
define 6^ ^4 = ^ J^i% dxa\{x)e-^^-\ The length 
scale L can be associated with the spatial extension of 
the initial pulse, see Fig. [l] The mean number of the 
wave packets in the state \a)r is given by TV = |ap. It is 
also convenient to introduce the dimensionless coupling 
strength 7 = irg^L and detuning 5 = {ko — A)L. To fa- 
cilitate computations we represent the initial state in the 
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where the displacement operators I?e,o are defined with 
help of mutually commuting operators b], and bj^ j,^^ , 
respectively. 

The statistics of the initial field, defining a probabil- 
ity Pain) to find n photons in the mode fcp, is given by 
the Poissonian distribution Pa{n) — e^^^. Due to 
the presence of the photonic dispersion, photons can leak 
from the mode fco to other modes on both branches of the 
spectrum by virtue of scattering processes, what induces 
nontrivial correlations between photons. Statistical prop- 
erties of the incident beam are thus being changed. A 
fraction of photons is refiected, and their statistics is also 
of great interest. We propose a calculation of the FCS in 
both forward and backward scattering channels, which is 
exact and thereby nonperturhative in both g and N . 

Defining FCS. — Generally speaking, FCS can be de- 
fined as a generating function F{X) — Y^^=o e'^"p(n) 
associated with a probability distribution p{n) of photon 
numbers n in some specific photonic state. The func- 
tion F{\) (InF(A)) generates m-th order moments (cu- 
mulants) of the distribution p{n): One has simply to con- 
sider the TO-th derivative with respect to at A = 0. De- 
pending on a physical meaning of the distribution p{n) , 
one can distinguish between different types of FCS. 



It is convenient to introduce the variable z 



anal- 



ogous to fugacity. In these terms, the Fourier expansion 
of the 27r-periodic function F{X) acquires a form of power 
series in 2, F(2) = X]^o •^"^'("') • As a function of z, F(z) 
can be analytically continued inside a circle of the unit 
radius \z\ < 1. In particular, one can treat z as a real- 
valued parameter on the interval —I < z < 1. Under this 
constraint the function F(z) is real-valued. The prob- 
ability normalization X]^o?'(") ~ ^(■^ = 0) = F{z ~ 
1) = 1 follows from the normalization of a quantum state 
whose statistical properties are being studied. 

In the FCS framework we want to characterize the 
scattering state |out) = 5*1 in). Due to the factorization 
property pi) this state can be represented as |out) = 
|out, a/'\/2% ® \a/V^)o, where |out,a/\/2)e is the scat- 
tering state in the model ([ij with the single (even) branch 
of the spectrum [231 US]; and \a/^/2)o is the coherent 
state in the odd sector. Note that in both cases the co- 
herence parameter equals a/\/2. 

While focusing mainly on the coherent initial state we 
define the momentum-resolved FCS of photons in the 
scattering state \out) corresponding to the initial state 

\in) = 



F^{Xr, A/) = (out|e' 



.iKNr.+iX,Ni 



out), 



(3) 
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with help of the wave packet number operators 



ji=0 



|0), ^ = rj (4) 



which possess nonnegative integer eigenvalues. Using the 
identity e''^f^« = 1 + ~ where — 

e^'^s, we express ([3| for |?/;) as 



E 



oo 
m=0 



(zr-lK(m)|so™p 



(5) 



+ ^ (z," - l)(zr - l)Pa(nK(m)|s„™|2 

n.m—O 

where the coefficients Snm quantify how much the FCS 
([3| differs from the initial, Poissonian one (for which all 
I Snm I = !)• The forward and backward scatterings are 
characterized by coefficients :— s„o a^id sjj := Som 
respectively. A computation of the coefficients exploits 
the explicit form of the |out) state ^23 . As a result, we 
find [22] the following expressions 
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where p = {S + i^)/2, Cn{p) = pe'^iju-iip) ~ Un(p)] and 
jn{p) are the spherical Bessel functions of the first kind. 
In the limits 7—^0 and 7 — >■ cx) we obtain sj^ = 1, sjj = 
(the whole pulse is transmitted) and sj'j ~ 0, s'„ = (—1)" 
(the whole pulse is reflected), respectively. In both cases 
the Poissonian statistics is retained. 
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FIG. 2: The probability distribution Pr(?^) for the number of 
particles in the mode fco in the forward-scattering channel, 
plotted for various values of coherent state parameter TV and 
the field-emitter interaction strength 7. Distributions for a 
different values of N are plotted in the Supplement. 




FIG. 3: The probability distribution pi{n) for the number of 
particles in the mode kg in the backward-scattered channel, 
plotted for various values of coherent state parameter N and 
the field-emitter interaction strength 7. Distributions for a 
different values of N are plotted in the Supplement. 



Results and discussion. — In Figs. [2] and [3] we plot the 
probabilities Pr,i{n) — Pct(n)|sJ^''p associated with ^ 
which describe the distribution of transmitted/reflected 
photons in the mode (—fco) after scattering. The mean 
numbers of photons corresponding to these distributions 
are presented in the Supplementary Material [55]. 

Analyzing the FCS we elucidate several important re- 
sults coming from the many-body character of scattering 
off the two-level system. First, the resulting distributions 
of photon numbers becomes essentially non-Poissonian 
when the interaction is switched on. Second, although 
the absolute value of density of particles in the mode fcg 
is reduced as a function of 7 this evolution is very non- 
trivial: PR{nko) is non-monotonous as a function of 7 
with a tendency towards forming a peak at fewer parti- 
cles number n while at some "f{N) a new peak emerges 
at some n* > N with a larger height. Third, the behav- 
ior of probabilities is weakly oscillatory as a function of 
n. One can understand these effects as a competition be- 
tween scattering into other-fc 7^ fco modes, backscattering 
and the initial Poisson distribution. 

To gain further insight we first analyze the situation 
of continuous laser radiation, L — > 00. This corresponds 
to p — ^ cx) with the constant ratio S/j. In this limit a 
number of universal results emerge which can be char- 
acterized in terms of transmission T = |tp and reflec- 
tion R = 1 — T probabilities. Here t — 6/{6 + i^). 
Analytical results for FCS (see Supplement for details) 
F{Xr,Xi) = 1 + e-^'fr + e-^'-fi + Ufi where - 
exp(A^5(e*^e - 1)) - 1 and N^=r,i = {NT,NR} suggest 
that the distribution for transmitted photons (A; — 0) 



is himodal, p{n) = (1 — e O'^n.o + e Pa,.{n) where 



lad 



N^. Similar results hold for reflected photons. 



Moreover, in this limit we are able to compute the FCS 
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for an arbitrary initial state. In particular, for the Fock 
(number) initial state \N) we again obtain a bimodal dis- 
tribution pAr(n) = Sn,o{^^T^)+Sn^NT^ in transmission 
channel. This result is highly non-classical and reflects 
the quantum nature of the number state: either all N 
photons are simultaneously transmitted or all reflected. 
Similar bimodal structure of the distribution takes place 
for the squeezed state. Exphcit result for the i^^(Ar,A;) 
for the arbitrary initial state \ip) is given in Supplement. 
The FCS in the right/left channels has an interesting fac- 
torization structure: the factors p,p{n) specific for the ini- 
tial state lip) multiplies by the factors |s„mp coming from 
the scattering itself. This leads to the bimodal structure 
of the resulting p{n). 

In the pulsed laser regime (finite L) the factors s„ have 
a quasi-periodic oscillatory (in n) behavior which can be 
approximated as sj^ ~ cos{^/yn/2) in the limit of large 
n. For generic n these functions of 7 and 6 are plotted in 
the Supplement. From this we infer that when the max- 
imum of the Poisson distribution Pa(n) (^ N) coincide 
with one of the maximum of the jsj^p function we get an 
enhancement of the scattering in the n* -particle channel. 
This manifestation of the many-body scattering effect is 
visible in Figs. ([2]) as a re-entrant effect described above. 
This phenomenon suggests an experimental way to tune 
to a resonance with n-particle scattering. We note how- 
ever that the available phase space volume for the irre- 
ducible n-particle scattering gets smaller as n increases 
(see e.g. }24j for generic arguments) and therefore many- 
body effects are suppressed by these natural causes. 

Scattering to the other fe-modes is not directly visi- 
ble from the FCS we computed here. In a sense, other- 
then-fco channels act as dissipation with respect to the 
ko channel. Careful analysis of the n-particle S'-matrix 
P5] reveals that a single-particle scattering to the fc ^ fco 
modes is suppressed by the energy conservation while es- 
sentially many-body scattering effects coming from the 
irreducible part of the n-particle scattering are generi- 
cally suppressed as (fc — fco)^" (this power is stronger 
at the resonance). This shows that while single-particle 
scattering to fc 7^ fcg is suppressed by energy conserva- 
tion, irreducible n-particle scattering is suppressed by the 
power of n and phase-space arguments. 

We comment further on the previous studies of statis- 
tics in this quantum optical setup. Indeed, after the work 
of MoUow there were many studies of phenomenon 
now referred as a resonance fluorescence. Statistics of 
classical light was studied long ago in [26] in a system of 
two-level atom coupled to a sm^te-mode photonic field. 
In this approach a strong coherent field approximation 
has been used (which results in replacement of operators 
by their coherent-state expectation values) as well as a 
Markovian approximation |27j . Here we avoid these ap- 
proximations. Moreover we can exactly evaluate several 
other FCS-like functions which clearly show generically 
non-Markovian character of a scattering in our setup. 



This studies will be reported elsewhere [25] . 

Experimental scheme. To measure correlations en- 
coded into F^{Xr, Xi) one needs, first of all, a frequency- 
resolved detection scheme. Second, we suggest to mea- 
sure the FCS using the philosophy employed in |16j in 
the cold atom context. Instead of studying the moments 
of distribution to a certain order one can directly plot 
histograms of measurement outcomes using homodyne 
detection scheme. Recent experiments on tomography 
of quantum states in multi-mode transmission lines cou- 
pled to q-bits have employed similar method [29| in the 
microwave domain. Typical interaction strength can be 
varied between 50 MHz -300MHz while the typical trans- 
mission line frequencies are of the order of 500MHz and 
environmental losses are typically small. We note that 
changing the pulse's duration one can effectively change 
the coupling strength 7. The line width of the source 
employed in |29j for microwave transmission lines cor- 
responds to the very large pulse duration limit (large 
L) where the universal effects described above can be 
directly observed. A selective enhancement of the n- 
particle scattering effects can be achieved via balancing 
between parameters of the initial state and field-emitter 
coupling strength. 

Conclusion. In conclusion, we have exactly computed 
for the first time a full counting statistics of a funda- 
mental problem in quantum optics - a pulse of coherent 
light propagating in a multi-mode channel and interact- 
ing with a two-level system. These results allow for the 
quantitative determination of the effects of many-body 
correlations coming from the interplay of photons' dis- 
persion, interaction of photons with emitter and physical 
parameters of a pulse (duration and average number of 
photons) . By tuning the parameters of a pulse one can se- 
lectively access n-particle correlations. For generic initial 
state the limiting form of statistics is given by a bimodal 
distribution. Effects described here can be directly ob- 
served experimentally with existing setups. 
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Supplementary Material 

Here we outline essential technical details of our derivations and present some additional results. 



THEORETICAL MODEL AND SCATTERING FORMALISM 



Theoretical model 



Our model is motivated by current experiments in nanophotonics where a one-dimensional (ID) transmission line 
is coupled to an emitter (an artificial atom, a quantum dot, a NV center) or cold atoms are coupled to evanescent 
modes. More generally, the model in question can be considered as a basic model of quantum optics, since quasi-lD 
models effectively emerge when the scattering is restricted to the s-wave channel. 

Its derivation relies on a number of approximations which are customary in quantum optics: (i) a dipole approxima- 
tion, (ii) the rotating wave approximation (RWA), and (iii) Born-Markov approximation. In addition, we linearize the 
photonic spectrum around some appropriately chosen working frequency flo which is commensurate with the emitter's 
transition frequency f2. Under these assumptions [except for (ii)] we obtain an effective low-energy Hamiltonian 



H 



E 



+ .90 ^ J dk{a'^f^ + a^fc)(cr"'" 



(SI) 



featuring the two-branch linear dispersion with right- = r = +) and left- = I = ~) propagating modes. To 
implement the RWA in a systematic way, we first perform the gauge transformation H — >■ WHU + i{dW / dt)U with 



U = exp 



which leads us the Hamiltonian 



H = ^ J dkS,ka^^f^a^k + ^(^^ + go ^ j dk{a\,,(j + a^fecr^ 
I dk{al,a+e^^"°'+a^ka-e-'^''°'), 
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(S2) 

(S3) 
(S4) 



where A = — J7o. As soon as 5o/^o ^ 1- the oscillating terms (S4| can be treated as a perturbation. In zeroth 
order they are simply neglected, what is equivalent to the RWA. 



Definition of wave packet field operators 

In order to define a coherent state we need field operators annihilating/creating states normalized to unity. The 
field operators Uk and a\, which are the Fourier transforms of a{x) and a^(x), do not suit for this purpose, as they 
fulfill the commutation relation [afc,a|,,] = 6{k — k'). To circumvent this difficulty, we construct wave packet field 
operators 

6fc = — / dxa{x)er'^'' , fot ^ — / dxa\x)e'^'' , (S5) 

which do satisfy the desired commutation relation [^/coi^Iq] — 1- They annihilate/create photons spatially localized 
on a finite interval of the length L. With their help we are able to introduce the coherent state 

t °° it" (h^ "l" 

=e"''»'''-o-"^o'"°°|0) =e-l"'°olV2y k,) ^gg^ 

^ — ' n\ 

n=0 



possessing the proper normalization {ako\otko) — 1. 
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Transformation to the "even-odd" basis 



Due to energy independence of the coupling constant go (which is the part of Born-Markov approximation), one 
can decouple the model defined by (S3) into two sectors. To this end, one introduces even (or symmetric) and odd 
(or antisymmetric) combinations of fields corresponding to the same energy 



_ flrfc + 0,1, -k _ Ork — ai,-fc 

Oek ! '^ofc 



V2 
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(S7) 



By virtue of this canonical transformation the Hamiltonian (S3) turns into a sum of the two terms, H = + iJ^, 
defined by 



= dk 



k al^Gek + 9 {a\f,cF + a^kCf^ 
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dkkal^,aok, 



(S8) 
(S9) 



where g = goV^- Note that the odd Hamiltonian Ho is noninteracting, and therefore odd modes do not scatter off a 
local emitter {So = !)• The even Hamiltonian H^, can be interpreted in terms of a chiral model with a single branch 
of the linear dispersion. Scattering in chiral models has been studied Ref. [53] for arbitrary initial states including 



the coherent state. If an initial state of the model (S3) defined in the original right-left basis admits a decomposition 



into a product state in the even-odd basis, and this is indeed the case for the coherent state 



(SIO) 



then it suffices to apply the results of [53] in order to find the scattering state Se\a/\/2)f,- All expressions necessary 
for this purpose are presented below. 



Scattering of the coherent state in the chiral model 



For the chiral model the outgoing state resulting from scattering of the initially coherent state \ak(,) ch (in the mode 
ko) has been constructed in Eq. (134) of Ref. ^23]. For completeness, we quote it here as well 



\ovLt,ako)ch = '5'ch|afco)ch, 

oo „ 



(Sll) 



e(L/2 > Xn> ■■■> X2> xi> -L/2) Y[ (l 



J(ko-Ag)Sxj + 



n 

+ e{L/2 >Xn> ...>X2> -L/2 > a;^)e-»(feo-'^9)(:^i+i^/2) (l _ ^riko-A,){x2+L/2)^ Y[ ' e'^^"~'^^'>^''^ + ''^ 



i=2 



where Sxj = Xj — xj^i > 0, Xn+i = L/2, as well as Ag = A — ing'^ and A = —2Trig'^ako/[{ko — Ag)\/L 
Let us now evaluate the factor 

T(vl^) = (O|e''^o^^o|out,afc„)ch = e-l"'=ol'/2+-Jo"fco (0|e''^o''^o ^^j^e-"^o^'=o |o) 



(S12) 



(S13) 



for arbitrary w^^^. Permuting e"'''o^'°o in (S13) to the right produces the shift of fields a^{xi) -> a^{xi) + ''^o^* 

spanning the scattering matrix 5ch. Reorganizing the resulting series for S'[w^^] = e^'=o'''=o S'che "'°o'''°o , we obtain the 
following expression 



SK] = S^[vl;L/2, -L/2] + S'[vl;L/2, -L/2]-^AI 



(S14) 



S2 



where 



and 



-L/2 



(S15) 



oo „ 

S''{vl^;x,y\ = d^ix-y) + ^" / rfa;„ . . . da;ie(x > x„ > . . . > a;i > y > -L/2) 

n=l 

X d^{x- x,,)a^x„)e'''°''"d^{xn - Xn-i)a\xn-i)e''"°''"-' ■ ■ ■ a\xi)e^'""'' d^{xi - y), 

OO „ 

S''[vl^;x,y] = dy{x -y) + ^y'' / dx„ . ..dxiQix > x„ > . . . > xi > y > -L/2) 

n=l 

X d^{x- x„)a^(x„)e'''''^"d„(2:„ - x,,^i)a\x„_i)e'''°''"-' . . . a\xi)e''""'^d^{xi - y). 



The operators (S16) and (S17) are defined with help of the kernels 

dv{x) = 
dv{x) = 



P+-P- 

(^-ip+x _|_ P+ ^-ip-x 
P+ -P- P+ -P- 



where p± = p±(w^g) are the roots of the quadratic equation 

+ {ko - Ag)p - 2TTg'^ak„vljL = 0. 
From this representation we straightforwardly get J-{v1^^) = dy{L) = dy{L, ak^). 



(S16) 

(S17) 

(S18) 
(S19) 

(S20) 



COMPUTING FCS 



Due to factorization of the initial state (|S10|) and of the scattering matrix S = SeSo = Se^ in the even-odd basis, 

(S21) 



we also obtain the scattering state in the product form 

|out) = 5|in) = |out,a/V2)e ® \a/V2)o 
As we want to measure moments of the operators — r,l), we start our consideration of the FCS from 



oo 
11=0 

n=0 \ / 



be.kp + £,bo,k o \ " 

V2 



(S22) 



It follows that 



iXrNr+iXiNi 



C30 OO 

n— m— 

oo 

+ E —,{e^''--l){e'''"'-l)brb]"'\0){0\bTb^ 



(S23) 



n,m— 



Representing — (6e,A;o + ^^omM^ calculating the quantum average of (S23t in the outgoing state (S211, 
we obatin 

oo ^ oo ^ 



n=0 



m=0 



oo 



nm I 1 



(S24) 



n.m— 



S3 



where 



9" 9" 



-|a|V2 



5(t;*)" 






Qm 






Qn 









|exp 



^/2 

^/2 , 



V2 , 



lout, aNI 



■u:=o 



(S25) 



Introducing new variables 'w± = f (i'^ ^D' obtain 

-|aP/2 ^a^"+™ { d _ d Y { d d 



Sn 



where 



and 



Then 



dw-^ dw- J \ dw-\- dw- 



o \ m 

—) I 



w±—Q 



1 / a 



2"+"' \dw 



71 / \ m 

a 

dw 



dw{L) 



, w = w+, 



a"+"s„™, (S26) 



(S27) 



io=0 



/ r /r>x f ■ X sin(KL/2) 

cos{kL/2) - i{5 + 17) — ^ — 

kL 



, nL = ^{6 + i-iY + 87U;. 



(S28) 



F„(A,,Az) = l + ^(e^^'-"-lK(n)|s„o|'+ ^(e^^""-lK(m)|so,n|' 

n— m— 

00 

+ el"l' ^ (e'^'-"-l)(e'^'™-l)p„(n)p„(m)|wr 



(S29) 



Introducing notations p = {5 + i7)/2 and t = —^w/p, and using an expression for the generating function of the 
spherical Bessel functions 



- cos(V'/32 - 2pt) = V -, j„_i(p), 



(S30) 



n=0 



we cast (S28 1 to 



dw{L) = Y] —,Cn{p), Cn{p) pe'''[j„_i(p) - ijnip)]- 



(S31) 



n=0 



and from Eq. (S27) explicitly calculate 



p=0 



Cp+m(p)- 



(S32) 



Note that these coefficients do not depend on a. For illustration we plot the coefficients s^^ and s„ in Figs. S4 



and 



At large |p| 3> 1 we use an asymptotic expression for c„ w i" and observe a factorization of the coefficients 
Snm = i^r™, where t = d/{S + 17) and r = —i"//(6 + ij) are transmission and reflection amplitudes of the single 
photon scattering, respectively. 



S4 
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FIG. S4: The dependence of magnitudes of coefficients sj^'', introduced at Eq. (|6| on n at various coupling constants 7 and 
detunings S. 




10.0 



FIG. S5: The dependence of magnitudes of coefficients s^' , introduced at Eq. Q on n at various coupfing constants 7, while 
the detuning 5 = 0. This is another representation of same data from Fig. |S4| 



S5 

Then, in the hmit L ^ oo (continuous laser radiation) we obtain 

oo oo 

F„(A,, AO = 1 + ^(e*^--" - l)p„(n)r" + ^ (e^^'" - lK(m)i?'" 

n— m— 

oo 

^ (e'^-" - l)(e'^'" - l)p„(n)p„(m)r"i?™. (S33) 



e 

n.m— 



After a simple transformation we arrive at the explicit formula 



n— m— 

^-^ n\m\ 

n.m— 



where T = \t\^ and i? = jrp = 1 - T. 

Considering the case A/ = we obtain the bimodal distribution for transmitted photons 

F(A,, 0) = 1 ~ e-^^ + ^-R^^TN{e^-^-i) ^ 
which interpolates between Poissonian and zero occupancy distributions. Analogously for reflected photons (A^ — 0). 



In the continuous laser radiation limit we can generalize (S33l for arbitrary initial state IV'fco) ~ S^o^"l")'=o 
created in the mode fco with the distribution of photon numbers v^n) = IV'nP- This becomes possible due the 
existence of the identity resolution 



1 = / — (S36) 

TT 



Then, 

— e-l"l'/2a"+'"(a|Vfco) = s„™ / — e-l"l'a"+™ V ^t^, = {n ^ m)\^^^^,,,^. (S37) 

and 

oo oo 

F^(A,, AO = 1 + y (e^^"-" - 1) p(n)r" + ^ (e'^'™ - l) p(m)i?™ 

n— m— 

+ £ ('''^"'Ve*^'-"-!) (e'^'™-l)p(n + m)r"i?", (S38) 

Thus, for the Fock state with iV photons = <5„.Ar we obtain 

Fjv(A„ 0) = 1 + (e^^-^ - 1) T^, (S39) 

which means that either all iV photons are coherently transmitted with the probability or all of them are reflected 
with the probability 1 - T^. 

For the squeezed state |C) = 5'(C)|0) where S'(C) = cxp(iC*^ls - K(^i?)^) ^"^^ ^ ^ l^l*^*^ coefficients V2n = 
(2n - l)!!(-e^^ tanh |C|)"/\/(2n)! cosh |C|, while V2«+i = 0. Then 

Fc(A„0) - (1 - rf^),5„,o + rfcPc'N (S40) 
where = cosh |C'|/cosh and p^in) is obtained from pQ{n) by replacing |C| by = arctanh(r^ tanh |C|). 



MOMENTS OF THE DISTRIBUTION 



Here we present some additional results on the mean number of photons remaining in the mode in the transmission 
and reflection channels as well as probability distribution function for different set of parameters. 
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FIG. S8: The probability distribution pi{n) for the number of particles in the mode ko in the /-channel, plotted for various 
values of iV and the interaction strength 7. 



